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Collective Mode in a (i^2_„2 + idxy Superconductor 



O 

o 

o 

(N 

o 



o 
o 

u 

Oh' 



X3 

o 

> 

m 

o 

On 
-(— > 



I 

O 

o 



X 



A. V. Balatsky^, P. Kumar^ and J. R. SchriefFer'^ 

^Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA. 

^ Department of Physics, PO Box 118440, University of Florida, Gainesville FL 32611 

^ NHMFL, Florida State University, 1800 E. Paul Dirac Dr. Tallahassee, FL 32310 

(Received February 1, 2008) 

We consider a superconducting state with a mixed symmetry order parameter components, e.g. 
d + is or d + id! with d' = dxy. We argue for the existence of the new orbital magnetization 
mode which corresponds to the osciUations of relative phase 4> between two components around an 
equilibrium value oi (f> = ^. It is similar to the so called "clapping" mode in superfluid ^He — A. We 
estimate the frequency of this mode luq {B, T) depending on the field and temperature for the specific 
case of magnetic field induced d' = dxy state. We find that this mode is tunable with an applied 
magnetic field with uio{B,T) oc BAo, where Ao is the magnitude of the d-wave order parameter. 
We also estimate the velocity s{B,T) of this mode. Pacs Numbers: 74.20.-z, 74.25.Nf 



The order parameter in high-Tc superconductors is 
widely believed K-M to be of a (1^2 _y2 symmetry. How- 
ever more careful consideration indicates that the sym- 
metry of the state in high-Tc might be lower in a num- 
ber of cases. The symmetry allows the secondary com- 
ponents to be generated whenever there is a perturbing 
field. These secondary components {is or id^y) genera- 
tion has been addressed in recent literature for the case of 
inhomogeneity due to wall scattering Iq-U^] or due to vor- 
tex texture ||l2|-[l6|. Similarly, the generated id' compo- 
nent near magnetic impurity and spontaneous violation 
of time reversal symmetry with global d + id' has also 
been discussed [[l7|-]l9[. Another possibility for a global 
d + id' phase has been pointed out po|-p2| where the ex- 
ternal magnetic field applied to two-dimensional d-wave 
superconductor generates id' . 

Direct experimental determination of the second super- 
conducting component nevertheless remains a substantial 
challenge, see e.g. g. The smallness of induced com- 
ponent and the large superconducting response of the 
underlying d component makes the detection of the sec- 
ondary component difficult. 

We present here an excitation which constitutes a di- 
rect test of the induced component of the order param- 
eter. This excitation is uniquely tied to the existence of 
the secondary out of phase component of the order pa- 
rameter. Consider the most general situation of a com- 
plex order parameter which can be generally written as 
Ao-|-iAi where Aq is the original dx2_y2 component and 
Ai is the induced s or dxy component which is orthogonal 
to the initial dx2-y2 state. Define the global (Josephson) 
phase of the order parameter u and the relative phase 
as: 



|Ao| + exp{i(j){v))\l^i\]exp{iv{Y)) 



(1) 



If the order parameter A is to be defined at the Fermi 
surface, then the functions A^, i = 0, 1 implicitly contain 
the angular dependence. The global phase v{y) can be 
position dependent and even singular as is the case for the 
vortex configuration. We will focus on the relative phase 
(j). Rs dynamics by definition is related to the appearance 
of the secondary component Ai. While the equilibrium 
value of the relative phase in most cases is argued to be 
±-| resulting in a violation of the time reversal symme- 
try, (f) may oscillate around this value and represent a 
new collective mode that can be used to detect the in- 
duced component. Similar relative phase mode has been 
considered before in context of two band superconduc- 
tors ll^ and heavy fermions [gj| and in [g5|, where the 
collective modes of d-wave superconductor in zero field 
were considered. 

This collective mode is akin to a spin wave in the or- 
bital moment of the superconducting ground state. The 
dispersion of this mode has a gap similar to the Larmour 
frequency in case of spins. As has been shown in [ p4[ , 
the relative phase oscillations also can be viewed as an 
internal Josephson effect, similar to the clapping mode 
in '^He — A p^ or to the relative phase dynamics in a 
2-band supercondutor |g^ . In the sense of a magnetic ex- 
citation, that can respond to a time dependent magnetic 
field in a resonant manner, this mode is also comparable 
to the longitudinal NMR in ^He - A [|6). The details of 
the dynamics are clearly model dependent. 

To be specific we will focus on the field induced Ai = 
idxy secondary component in the bulk. In this case quasi- 
particle spectrum is fully gapped in the field with the 
minimal gap vanishing at zero field. The d -(- id' state 
breakes time reversal symmetry and has a finite mag- 
netic moment M^ perpendicular to the superconducting 
plane p]],p3 . The mode discussed above in this case will 
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correspond to the longitudinal oscillations of the conden- 
sate magnetic moment around its equilibrium value. The 
frequency of this mode will contain information about 
the details of the coupling between components and will 
provide insight into the overall order parameter struc- 
ture. Below we will focus on a 2d superconductor at 
the fields H ~ Hc2 (T) and therefore we will ignore small 
{0{H—Hc2)) difference between induction B and applied 
field H||z. 

We find the resonance frequency for the relative phase 
oscillations of (f>. It turns out to be a gapped propagating 
mode with: 



We now will prove the existence of the "clapping mode" 
for d-\-idl state. The free energy has the standard form : 



F = F. 



^o = |^|AoP 
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uj^{B,k) = ujI{B,T) + sl{B,T)hkj, 



MB) 
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W) 



MB,T) 
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(2) 



(3) 



Here s = {a + bB^)AQ is the mode velocity, a,b are 
some constants and ry is a constant, a measure of the 
strength of the interaction, which we discuss in more 
detail below; 7V(0) is the Density of States (DOS) at 
the Fermi level; both i,j refer to in-plane coordinates 
x,y. This mode is a longitudinal magnetization oscilla- 
tion 6Mz{t, r) DC dMzexp{iuj{B, k)t — ikr) and is tunable 
by the external field. There are, in effect, two conse- 
quences arising from a non-zero rj which may be used to 
estimate its magnitude. In the presence of a secondary 
order parameter Ai = id^y, the temperature dependent 
upper critical field Hc2{T) has an additional contribu- 
tion which is quadratic in (1 — T /T^). This results in an 
upward curvature with a scaling field which depends on 
ry. Thus the upward curvature in Hc2{T) and the reso- 
nance frequency should be related in that they are both 
derivable form the same energy scale. 






clapping mode 
similar to 3He-A 
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where Aq, Ai are d and d' components of the order pa- 
rameter, oo — ao(T/ro — 1) , To is the ordering temper- 
ature for Aq. The corresponding ai > = iVo is always 
positive, as are 6o > and Kij,Klj > 0. Di — di — i^Ai, 
with B = V X A. For simplicity we take Kij — KSij and 
similarly for K' term. Because we are near Hc2, the 
flux expulsion is nearly non-existent and B and H are 
nearly equal. The interaction term in this form has been 
proposed earlier |20[|2^ . One can see that the GL free 
energy is indeed a scalar. The interaction term, using 
[Dx, Dy] — ieB can be written as 



-^7? 



leBJAoAl 



h.c. 



(8) 



which corresponds to a coupling of the magnetic field 
with an intrinsic orbital moment along z-axis:< M^ >— 
qAoAl + h.c. 

Up to now the work on secondary component has been 
focused on the equilibrium solution for a particular model 
for induced component, say in d-\-is or d+id' state. Here 
we will assume that the amplitude for the secondary Ai 
has been developed and we will look at the relative phase 
(j) oscillations only. It is convenient to introduce the re- 
spective phases of each component: 



Ao = \Ao\exp{i(f>o), Ai = |Ai|ea;p(z0i) 



(9) 
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FIG. 1. Schematic representation of the longitudinal mag- 
netization oscillations caused by the oscillation of the relative 
phase angle between Ao and Ai with an equilibrium value 
7r/2. The frequency of this mode is linear in the B field and 
is also linear in the magnitude of the Ao . This mode is there- 
fore tunable by the field and by temperature. We estimate 
ujo{B,T) ~ Ai(_B,r) and is in the Kelvin range. 



which are related to the introduced above v = (j)o,4> — 
(/>i — (/iQ. Derivation of the Eq.(g) then proceeds as follows. 
We use the Josephson relations for each component: 

OF ■ OF 



ON, 



fJ-i, Ni 



-,i^O,l 



(10) 



Where fii are chemical potentials for particles in Aq, Ai 
condensate and similarly, Ni are conjugated number of 
particles. In the equilibrium, where /io = /^i = A*, wc find 
for a relative phase motion: 



r dfiQ dfii _ dfiQ dF 
^dNo dNi dNihcl) 



(11) 
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where in the above we have taken into account the 
fact that although jjq = /xi their derivatives are differ- 
ent. The term in the brackets is on the order of iV(0): 
[|^ + l^-2|^] - P~' - ^(0)- For the general values 
of the relative phase </> the terms in free energy F that 
contribute to the derivative in Eq.(pl|) are Fint, Eq.(|^) 

K' 



and terms with gradients Kij , _ 

Fvnt = 77-62 1 Ao 1 1 All sin ( 



(12) 



with minimum F reached at (f) = —t: /2sgn{Bz). In 
general, there is an additional term proportional to the 
squares of the two (the bare and the induced) order pa- 
rameters. However it contributes little new to the physics 
of the problem. Its effects are largely quantitative. We 
stress here that the phase (/)o of Ag is not constant in the 
external field in the mixed state. Nevertheless the F^^t 
in Eq.(n2) depends on the relative phase only. 




FIG. 2. The Fint profile as a function of the relative phase 
angle between Ao and Ai is shown. Although the minimum 
is reached ai tf) — ^ the finite stiffness for (f> oscillations leads 
to the finite frequency mode at ujo{B,T). We have ignored 
the dependence on the sign of B in the figure. 



We find: 



-TylAollAilB^cos^-s^V^ 



(13) 



Here s is given by Eq.(P). Minimizing the free energy 
Eq. (Hjpy^) with respect to the magnitude from this equa- 
tion, Eq.(p|) follows immediately. The velocity s is field 
dependent and also can be used in experiments to detect 
the presence of id' component. 

We can estimate the magnitude of the energy gap up 
to a 0(1) prefactors: 



wo{B,T)~\Ai{B,T)\ 



rjB 



Ao(S,T) 



(14) 



as one can easily see from minimization of the total F 
with respect to Ai. For estimate for the coupling con- 
stant 77 and magnitude of Ai see [g2|. We estimate thus 



MB,T) 



(-)(—) 



(15) 



This result turns out to be very similar to the "clapping 
mode" in '^He — A, where frequency was found to scale 
with the gap in the whole temperature range as well PSJ . 
We also find asymptotics in H: 



UJo 



H/Hc2 H < Hc2 
\/Hc2 — H H -^ Hc2 



(16) 



and we note immediately that the mode frequency ujq <C 
Aq Therefore this mode will be sharp. The damping 
coming from the low-lying quasiparticles in the nodes of 
d-wave gap will not affect this mode because the phase 
space for the decay will be small. 

In above estimates we assumed that the field induced 
component id' will be present even outside the GL region 
for which the above formulas have been derived. We as- 
sumed that once the field induced component is present 
in GL region it will persist to a lower temperatures. Real- 
izing the limitations of estimates made beyond GL regime 
we present Eq.(nq) as an order of magnitude estimate 
only. From Eq.(|15|) it follows that the gap cjq can be 
made in the range of 0.1 - IK {2 - 20GHz) in the field 
H = 1- lOT similar to the Ai estimates ^. 

The longitudinal oscillations of magnetization M^ will 
also lead to the resonance at ujo{B,T) in the AC sus- 
ceptibility of the superconducting state. Experimentally 
the proposed clapping mode can be observed in NMR 
in the field. Another possible method of detecting this 
resonance is the ultrasound attenuation (taken at lower 
fields), used previously in heavy fermions to investigate 
the excitations of the order parameter. For any of the 
resonance techniques used to search for the "clapping 
mode" resonance at uJoiB,T) Eq.([l5|) the vortex cores 
contribution would be the biggest source of background. 
We note however that ujo{B,T) can still be measured if 
the resonance is sharp enough. 

In conclusion, we considered the relative phase oscilla- 
tion mode that can be used to detect the secondary com- 
ponent in the time reversal violating superconducting 
state, such as d + id' or d+is. Two components of the or- 
der parameter are characterized by respective amplitudes 
and phases. The relative phase between two components 
can oscillate around its equilibrium value = ±^. This 
mode is very similar to the "clapping mode" in super- 
fluid ^He — A. For a specific model we choose the case of 
field induced id' component. We show how the relative 
phase mode frequency is governed by external magnetic 
field and d-wave gap magnitude uJo{B, T) oc BAq and is 
therefore tunable by external magnetic field. The scale 
for 0^0 and mode velocity s is set by the magnitude of 
the induced gap Ai and we estimated the frequency for 
different fields, Eq.dlS). Among other probes this mode 
could be experimentally detected by investigating the ac 
magnetic susceptibility and possibly by ultrasound atten- 
uation in the in the mixed state as a function of applied 
field H. 
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